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We study the weakly disordered Bose-Hubbard model on a cubic lattice through a one-loop 
renormalization group analysis of the corresponding effective field theory which is explicitly derived 
by combining a strong-coupling expansion with a replica average over the disorder. The method is 
applied not only to generic uncorrelated on-site disorder but also to simultaneous hopping disorder 
correlated with the differences of adjacent disorder potentials. Such correlations are inherent in 
fine-grained optical speckle potentials used as a source of disorder in optical lattice experiments. As 
a result of strong coupling, the strength of the replica mixing disorder vertex, responsible for the 
emergence of a Bose glass, crucially depends on the chemical potential and the Hubbard repulsion 
and vanishes to leading order in the disorder at commensurate boson fillings. As a consequence, 
at such fillings a direct transition between the Mott-insulator and the superfluid in the presence of 
disorder cannot be excluded on the basis of a one-loop calculation. At incommensurate fillings, at 
a certain length scale, the Mott insulator will eventually become unstable towards the formation of 
a Boss glass. Phase diagrams as a function of the microscopic parameters are presented and the 
finite-size crossover between the Mott-insulating state and the Bose glass is analyzed. 



PACS numbers: 67.85.Hj, 03.75.Lm, 72.15.Rn 



I. INTRODUCTION 

The disordered Bose-Hubbard (BH) model can be 
viewed as the fruit fly for studying the interplay be- 
tween disorder and strong interaction physics. In the 
absence of disorder, two phases compete, both of which 
have a gap in the single-paricle spectrum and hence 
are incompressible: the Mott insulator (MI) which wins 
when the on-site repulsions dominate and a superfluid 
(SF) which obtains when the kinetic energy exceeds a 
critical threshold. Disorder introduces a new avenue 
for localization independent of the interaction strength. 
Fisher, et ali argued that a Bose glass (BG), a gap- 
less insulating state, always intervenes upon the disrup- 
tion of the SF in the presence of disorder. Nonethe- 
less, this finding remains controversial as simulations 
and analytical arguments both support 2 -^^^^ and 
negate^^^^^^^^^ this claim, the latter fo- 
cussing predominantly on the case of commensurate bo- 
son fillings. In the absence of disorder, the suppression of 
density fluctuations at commensurate fillings is known to 
change the nature of the transitions from mean-field like 
to the universality class of the (D + l)-dimensional XY 
model 20 ' 3 - with D the spatial dimension, suggesting that 
the system might respond differently to an infinitesimal 
amount of disorder. It has been argued by Fisher, et ali 
that such a scenario is in principle possible but unlikely 
since disorder is expected to destroy commensuration. 

Certainly, the experimental realization of the MI-SF 
transition in a gas of ultracold bosonic atoms in an opti- 
cal lattice 2 ^ and the possibility of introducing disorder 
in a highly controlled manner^ e.g. by superimpos- 
ing an optical speckle potential , 24 i 25 i 26 i 27 nas triggered 
hope that experiments might address the questions raised 
above. However, the experiments suffer from serious limi- 



tations, such as, inhomogeneous densities due to the pres- 
ence of an optical trap potential, heat generation by the 
lasers, and small system sizes, which might ultimately 
render them incapable of answering questions about uni- 
form systems in the thermodynamic limit. This raises the 
question as to whether or not optical lattice experiments 
will live up to the expectation that they provide the ulti- 
mate quantum simulators of paradigmatic Hamiltonians 
as the one of the disordered BH model, 
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This model describes bosons with corresponding creation 
(annihilation) operators b\ (5j) hopping with (disordered) 
amplitudes ty between nearest-neighbor bonds of 
a D dimensional hypercubic lattice in the presence of 
an on-site disorder potential ej and subject to a local 
Coulomb repulsion U. Further, n, = b\bi denotes the oc- 
cupation number operator and fj, the chemical potential. 
In the following we revisit the general disorder problem 
with an eye for the special form of disorder that is in- 
duced by a fine-grained optical speckle lens used in recent 
experiments^ 

Detailed measurements of the speckle potential 2 ^ and 
subsequent calculations based on the imaginary time evo- 
lution of localized Wannier states in the presence of the 
measured potentia l 27 ' 28 have revealed the induced dis- 
tributions of the microscopic parameters. The potential 
disorder has been shown to follow the asymmetric dis- 
tribution P(ei) = e- £ -/ A /A for a > and P(e») = for 
ej < 0. The reason for the one-sidedness of the distribu- 
tion is the use of a blue-detuned speckle potential which 
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can only increase the potential energy of a lattice site. 
Since correlations between different sites are small even 
for nearest neighbor sites, we will assume the potential 
disorder to be uncorrelated. Simultaneously, the speckle 
field induces disorder in the hopping amplitudes which is 
correlated with the difference of the disorder potentials 
on adjacent sites such that ty — t + 7(e, — £j) 2 - Here, 
t denotes the nearest-neigbor hopping in the absence of 
disorder. The relative width of the distribution of the 
on-site repulsion U has been shown to be negligible. 

In this work we augment a strong-coupling expan- 
sion around the localized limi t 29 i 30 with a replica aver- 
age over the disorder. In the clean limit, this procedure 
captures the quantum phase transition between the MI 
and the SF. The resulting effective field theory is then 
studied within a one-loop renormalization group analy- 
sis. The combination of strong-coupling expansion and 
replica calculus outlined in this work allows for a system- 
atic calculation of thermodynamic phase diagrams of the 
disordered BH model. 

Our major finding is that the effective replica-mixing 
disorder vertex, which is responsible for the emergence of 
a BG, depends crucially on the value of the chemical po- 
tential. Irrespective of the special form of the disorder, on 
approaching commensurate fillings, the disorder vertex is 
strongly suppressed and vanishes to leading order given 
by the variance of the disorder distribution. As a conse- 
quence, the length scale at which the MI becomes unsta- 
ble towards the formation of a BG becomes increasingly 
large. On the basis of a one-loop renormalization-group 
analysis which consistently treats the leading disorder 
contributions a direct MI-SF transition remains possi- 
ble in the presence of disorder at commensurate fillings. 
Whereas, as we show, a conclusive answer necessitates 
a considerably more complicated 2-loop calculation, the 
anomalous suppression of the generic disorder vertex and 
the presence of additional disorder operators which turn 
marginal at commensurate fillings strongly suggests that 
the nature of the transition at commensuration is in a 
different universality class from the transition at incom- 
mensurate fillings. 

The outline of the paper is as follows. In Sec. [H] we 
derive the effective field theory for the disordered BH 
model in the strong-coupling and weak disorder regime, 
providing explicit expressions relating the effective cou- 
pling constants to the microscopic parameters. To lay 
plain the details of the derivation, we start with a brief 
review of the Hubbard-Stratonovich transformation un- 
derlying the strong-coupling expansion in the clean limit 
(Sec. Ill A[) and illustrate the nontrivial coupling of on- 
site potential disorder to the dual theory. Employing 
the replica trick to average over the disorder, we ob- 
tain the effective field theory (Sec. Ill Bp . In the case 
of simultaneous hopping and on-site disorder, we derive 
the effective field theory by first performing a disorder 
average of the replicated system and by then general- 
izing the consecutive Hubbard-Stratonovich transforma- 
tion and strong-coupling expansion (Sec. Ill Cp . Starting 



from the effective action, we derive in Sec. Mil the one- 
loop renormalization-group equations. We present the 
integrated numerical solutions in Sec. IIII Bl which en- 
able a direct study of the instability of the MI towards 
the formation of a BG. The resulting phase diagrams are 
presented in Sec. IIVI Finally, in Sec. El we summarize 
our results, compare with other analytical work, and dis- 
cuss the relevance of our findings to recent optical lattice 
experiments. 



II. EFFECTIVE FIELD THEORY 

In this section, we derive the effective long-wavelength 
field theory for the weakly disordered BH model in the 
strong-coupling limit. We focus on temperatures much 
smaller than the on-site Hubbard repulsion, /3U <C I 
where (3 denotes the inverse temperature. In order to 
set up a field theoretical description, it is useful to ex- 
press the partition function of the disordered BH model 
([T]) as a path integral in imaginary time r G [0; 0) over 
complex bosonic coherent states, bi\4>i) = 4>i\4>i), 
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where we have explicitly built in the correlation between 
the on-site disorder potential and the hopping disor- 
der SUj. Note that in this model a perfect correlation 
between the on-site and hopping disorder is assumed. 
This seems reasonable for disorder induced by an optical 
speckle lens since all distributions of the microscopic pa- 
rameters as well as the correlations between them origi- 
nate from the same speckle disorder potential. A pos- 
sible generalization would by a model with two sepa- 
rate distributions -P(ei) and P(Stij) with tunable cross- 
correlations. 



A. Clean system 

We start with a brief review of the derivation of the ef- 
fective strong-coupling theory in the clean system (e, = 
on all lattice sites) . The formation of the MI state takes 
place in the regime where the Coulomb repulsion U dom- 
inates over the hopping amplitude t. Therefore, it is de- 
sirable to find a dual description in which residual vertex 
corrections are controlled by the smallness of t/U. Such 
an expansion around the localized limit can be achieved 
by a Hubbard-Stratonovich transformation as sketched in 
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Refs. 1 2!) and outlined in detail in Ref. [30]. The cru- 
cial step consists in the decoupling of the boson hopping 
term by introducing a complex auxiliary field ipi(r), 



= / drV (0*^ + 0^*), (3) 
Jo 



which is easily verified by completing the square and in- 
tegrating over the fields ipi(r). Here, T -1 denotes the 
inverse of the hopping matrix which is given by tij = t 
for i,j nearest neighbors and — otherwise. After 
taking the trace over the original bosonic fields & (r) , we 
rewrite the partition function of the BH model as 



Z = Z J V[^,i)*]e' [s '° )+S '^ 

5; = -ln(T T exp[/ dr£(^(T)St(T)+h.c.)]>o, (4) 
Jo 

where the average (. . .)o has to be taken with respect to 
the local on-site Hamiltonian 7io = ^i[^hi(hi — l)— iihi] 
which is diagonal in the local occupation number basis 
|{rii}). Further, Z Q denotes the partition function in the 

localized limit, Z Q = Tve- pn °. Note that we re-expressed 
the average over bosonic fields 4>i(r) as an operator av- 
erage where &J(t) = exp(7i i")S- exp(— Hot) and T T , the 
time-ordering operator. 

In principle, the part S'^ can be expanded to any de- 
sired order in the fields ip since the coefficients are simply 
related to bosonic Green functions of the local Hamilto- 
nian 7io- However, already the calculation of the coeffi- 
cients of the quartic terms ~ tp is tedious since it neces- 
sitates the evaluation of the two-particle bosonic Green 
function. Note that H.q represents an interacting local 
problem and consequently Wick's theorem does not ap- 
ply. After expanding SL up to quartic order, we take 
the continuum limit and perform a temporal and spatial 
gradient expansion to obtain the effective action 
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Here, the spatial gradient term with coefficient K% 
1/ {zt) results from the long- wavelength expansion of 



(0) 



and z — 2D denotes the coordination number of the D- 
dimensional hypercubic lattice. The mass coefficient of 



the theory has contributions from both S, 
is given by = K { 3 0) + C?(°) where 
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and S'^ and 
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is identical to the Fourier transform of the single parti- 
cle Green function Gi(r — r') = — (T T Si(r)S|(r'))o in the 
zero- frequency limit. In order to calculate G;(t), we first 
realize that at sufficiently low temperatures (j3U -C 1) 
the partition function is dominated by the ground-state 
configuration m — to where to is the integer minimiz- 
ing the site energy e„ = —/in + ^n(n — 1), hence the 
smallest integer larger than n/U. In this approximation 
the single-particle Green function is easily calculated by 
inserting a complete set of states, 1 = J2 n l n )( n l> 



Gi(r) = -(m+l)e- £ + r e(r) 
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Here Q(x) — 1 for x > and 0(x) = for x < de- 
notes the Heavyside function. For brevity, we have de- 
fined e + = e m+ i - e m = mU - /i and e_ = e m -\ - e m = 
(I — m)U + /1. Using this result, we obtain the mass 
coefficient in the zero-temperature limit, 
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The coefficients of the temporal gradient terms are given 
by derivatives of the mass coefficient with respect to the 
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chemical potential, k[ ^ = 

which directly follows from the temporal gradient expan- 
sion but can also be shown to be a consequence of the 
U(l) gauge symmetry. 29 

The interaction vertex is given by the connected parts 
of the two particle Green function in the static limit. 
In order to calculate the two particle Green function, 
GY{t u t 2 ,t 3 ,t 4 ) = (T T bi{n)bi{T2)bl{T 3 )bl(T 4 )), for every 
possible time ordering we simply insert identities 1 = 
En an< l evaluate the resulting products of matrix 

elements of bosonic creation and annihilation operators. 
The calculation is tedious but straightforward^ Here we 
provide only the final result which is given by 
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with e± as defined above, e +2 — £m+2 _c m = {l + 2m)U — 
2[i, and e_ 2 = e m -2 - e m = (3 — 2m)U + 

To summarize, we have derived an effective strong- 
coupling field theory providing the explicit dependence of 
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FIG. 1: Mean-field phase diagram of the clean BH model as 
a function of fi/U and t/U showing the first three Mott-lobes 
with m = 1,2,3 bosons per lattice site. The phase bound- 
aries are obtained by the sign change of the mass coefficient 
(/x, t, U). The dashed lines indicate the tip positions of the 
lobes at which the linear time derivative term in the effective 



field theory vanishes, K\ 
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the effective parameters on the microscopic ones, namely 
the on-site Coulomb repulsion U , the hopping t, and the 
chemical potential /i. Note that in this dual representa- 
tion, the vertex correction H 1 - ' is of relative order t/U 
as a consequence of the underlying expansion around the 
localized limit. We have neglected vertex corrections con- 
taining additional temporal derivatives. 

At mean-field, the phase boundaries between the MI 
states with fillings m and the SF are obtained by the 
sign change of the mass coefficient R\°\fj,,t, U). The SF 
corresponds to the broken-symmetry phase (R^ < 0) 
where the SF density is given by p s ~ (b) = K^°\tp}^ 
The MI states are obtained for R(°> > where the mass 
coefficient corresponds to the Mott gap. We show in 
Fig. [1] the mean-field phase diagram as a function of 
the dimensionless parameters \i/U and t/U. The tips of 
the different Mott lobes at (fj,/U) m — yjm(m + 1) — 1 
correspond to commensurate boson fillings which greatly 
stabilize the MI states. 

Note that the coefficients and are given by 
first and second derivatives of the mass coefficient 
with respect to the chemical potential and are therefore 
related to the slope and curvature of the phase bound- 
aries. At the tips, we find K± = and > yielding 
a field theory with dynamical exponent d z = 1, whereas 

at incommensurate fillings K~[ ^ ^ and d z — 2. From 
simple power counting, we find that at incommensurate 



fillings, the interaction vertex is an irrelevant perturba- 
tion for D + 2 > 4 where D denotes the spatial dimension 
and therefore, for D = 3 the transitions are mean field. 
However, despite being irrelevant, the interaction vertex 
leads to a renormalization of the mass coefficient and 
therefore to a slight shift of the phase boundaries as we 
will see in Sec. ITVl 



B. Pure on-site disorder 



As a next step, we use the results of Sec. Ill Al to de- 
rive the effective field theory in the case of weak on-site 
disorder with a probability distribution P{ti). With- 
out loss of generality we can assume that the distribu- 
tion has zero mean, e = J deP(e)e = and variance 
e 2 = J deP(e)e 2 = A 2 . In case of distributions with non- 
zero mean such as that induced by an optical speckle 
field, we simply shift the distribution and redefine the 
chemical potential as ~p — fi — e. In the following we fo- 
cus on weak disorder, A/U -C 1 and therefore neglect 
contributions from higher moments of the distribution. 
Further, we assume the disorder to be uncorrelated on 
different sites which is a good approximation in the case 
of fine-grained speckle disorder. 27 i 28 

Since disorder enters only in the form of random shifts 
of on-site energies, the Hubbard-Stratonovich transfor- 
mation outlined in the previous section Til Al can be per- 
formed in exactly the same way, yielding the action 



[ y i 



+ K 2l \d T ^\ 2 + Gil^f + H^l 4 } , (10) 



where the coefficients are given by the ones of the clean 
system but with the chemical potential shifted by the 
disorder potential on the corresponding site. The con- 
tribution to mass term takes the form 



g (Q) (p-^u) 



m + 1 



(I — m)U + fi — ei mU — /i + e. 



.(11) 



From the above equation it is clear that the coupling to 
disorder is radically different from a generic random mass 
problem where one assumes disorder to couple linearly to 
the mass of the theory. Due to the strong-coupling ex- 
pansion, the effective mass setting the Mott gap depends 
in a complicated way on the chemical potential of the BH 
model and therefore on the disorder. 

In order to restore translational symmetry and to per- 
form a continuum limit, we employ the standard replica 
trick to average over the disorder. To calculate the 



5 



disorder-averaged free energy, one has to average the log- 
arithm of the partition function, In Z = lim n _ > o(-2' n — 
l)/n. Therefore, we introduce n replicas of the system, 
ipi(i~) with a = 1, . . . , n, the replica index, and define the 
effective, disorder averaged action S e s as 



P[V^Cle~ Seff[W>r}1 , (12) 



where (• • • ) denotes the average over the disorder. The 
cumulant expansion through second order can be written 
compactly 

a ab 

in terms of S^ is = Sdis[i>i, ipf *]• After taking the contin- 
uum limit and performing the spatial gradient expansion, 
we find that 

S cS = /"dT f d^V^V^Va+i^lcW'al 2 
+X 3 |VVa| 2 + R\i>a\ 2 + 

+ G E/ fd D r\Mr)\ 2 \Mr')\ 2 , (14) 

ab 

where we have again neglected all vertex corrections con- 
taining imaginary-time derivatives. The effects of the 
potential disorder are twofold. First, the mass and inter- 
action vertices are renormalized as 



R 
H 



K. 
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- 17(0) 



1 d 2 H^ 
' 2 <V 



A 



(15) 
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respectively. As before, the coefficients of the temporal 
gradient terms are given by derivatives of the mass coef- 
ficient with respect to the chemical potential, K± = — ^ 

and i\2 = Since the disorder does not couple to 

the hopping-matrix elements, the spatial gradient term 

remains unchanged, K% — = l/(zt). Second, dis- 

order generates a vertex correction which mixes different 
replicas and is non-local in imaginary time, generic for 
effective disorder terms on the level of the replica aver- 
aged effective action. The strength of the disorder vertex 
is 



G 



K 



(0) 



A 



(17) 



As it should, in the clean limit (A = 0) the replica 
mixing term is absent and the effective action (|14p re- 
duces to n identical copies of that for the clean system. 



Note that the effective disorder strength G is not sim- 
ply given by the variance A 2 of the microscopic disorder 
but is also a function of the microscopic parameters \i 
and U entering via the factor k[°\ As we have seen 
above, this is a consequence of the underlying strong- 
coupling expansion and the resulting non trivial coupling 

to disorder. As we have realized before, k[ ^ 



<9/i 



has a simple geometrical meaning and is related to the 
slope of the phase boundaries between the MI and su- 
perfluid states. To be more precise, is the projec- 

tion of the gradient of R^ on the chemical potential 
axis. Since the gradient is perpendicular to the lines of 
constant mass, that is, the mean-field phase boundary 



0, K\ 



(0) 



changes sign at the tip posi- 
1 of the Mott lobes. 



given by R(°> 

tions (fi/U) m = ^Jm(m + 1) — 1 of the Mott lobes. As 
a consequence, the disorder vertex G vanishes at such 
commensurate fillings. 



C. Speckle disorder 

In this section, we proceed to derive the effective 
field theory for the BH model in the presence of weak 
speckle disorder. We assume a one-sided distribution 
P( €i ) = e~ ei / A /A for > and P(e, t ) = for a < 
of the disorder potentials as extracted for optical speckle 
disorder— and assume that the disorder potentials are 
uncorrelated on different sites. Since the speckle poten- 
tial not only introduces disorder in the on-site energies 
as studied in the previous section but also in the nearest- 
neighbor hopping amplitudes, Stij — 7(5^ — ej) 2 , it is not 
possible to decouple the boson hopping term before per- 
forming a disorder average. Such a procedure would re- 
quire the inversion of a disordered hopping matrix. Alter- 
natively, we can first perform a disorder average and then 
generalize the Hubbard-Stratonovich transformation to 
decouple terms which are not site-diagonal. 

Starting from the coherent state path-integral repre- 
sentation ^ of the speckle disordered BH model, we 
again employ the replica trick to obtain the disorder- 
averaged action 
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where we have used the disorder averages ejCj — A 2 (l + 
Sij) and (ei - e 3 ) 2 e k = 2A 3 (1 + 5 lk + S jk ) for i ^ j, 
which are easily derived from the moments e™ = n!A" 
since the potentials are uncorrelated on different sites. 
Further, we have defined the disorder-shifted chemical 
potential and hopping as /I = A 4 — &T = £t — A and t = 
t+7(ej — £j) 2 = A + 27A 2 , respectively. Here, 5^ results 
from the on-site disorder. Since to order A 2 the coupling 
of the disorder to the hopping term has only the trivial 
effect of shifting the nearest neighbor hopping amplitude, 
we have included the next order (A 3 ) replica mixing term 
which corresponds to nearest neighbor hopping coupled 
to the occupation on adjacent sites. 

Despite the more complicated structure due to this 
contribution, it is still possible to decouple the boson 

(2) 

hopping terms by introducing auxiliary Hubbard- 
Stratonovich fields ipf (r) as 



c( 2 ) / l c(°) i ct 1 ) i c( 2 h 

-oj. _ ; nr./.et „/,a* 



4° } = EM^V r ^' 

ija 

S S = -E/ dr(^^r+h.c), 

ia 

E/ (^ a (T)C(T)|^(r')| 2 + h. 
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a , 



where T denotes the inverse of the hopping matrix 
T which has non-zero elements t for nearest neighbors. 
After taking the trace over the boson fields 4>f(r), 4>°;*(t), 
we write the effective action formally as 



by i T (4 1 \5W) 4 )o-i(5«(5W) 2 )o((5W) 2 ) . There- 
fore, one has to insert two additional occupation number 
operators hf(T), tj|(t') into the operator products mak- 
ing the calculation of the average more tedious. Likewise, 

the terms of order A 3 contain exactly one factor 5^ 
requiring the insertion of one extra occupation number 
operator. Since the calculation is lengthy but straight- 
forward and not particularly insightful we will only give 
the final results. 

After taking the continuum limit and performing a 
temporal and spatial gradient expansion, we obtain an 
effective action identical to Eq. |TD] but with a mass coef- 
ficient 
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TO + 1 



(o)7A 3 



i^ 0) A 2 , 



(21) 



where R^°\ k[ q) , and are identical to the coeffi- 

cients of the clean system but with the hopping ampli- 
tude and chemical potential replaced by the disorder av- 
eraged values, t = t + 27 A 2 and ~p = n — A, respectively. 
As before, m labels the different Mott lobes, K\ 



K 2 = 



dR 
dfi ' 

If^r, and K 3 = l/{zt) with z = 2D the co- 
ordination number. The dependence of the coefficients 
on the chemical potential and the on-site Coulomb re- 
pulsion enters via the functions e_ = (1 — m)U + ~p and 
e + = mU — Jl. For the interaction vertex, we find 
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where the average 
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>o over the bosonic fields 4>f(r) 

has to be taken with respect to the action cor- 
responding to the site-diagonal replicated Hamiltonian 



W$ n) = £J-7«ft? + - 1)]. As for the 

clean system, the averages in the cumulant expansion 
are conveniently calculated by using the correspondence 
with time-ordered products of boson operators &| (t) — 

exp(7^ n) T)S;| exp(— T-Cqt) and bi a (r). These can easily 
be evaluated by inserting a complete set of bosonic states 
in the occupation number basis. The terms of order A 
are again related to the one and two-body bosonic Green 
functions and the calculation is completely analogous to 
the one outlined in section Hi Al 

Extra complications arise in the calculation of the dis- 
order terms. The contributions of order A 2 from the 
on-site disorder are given by averages containing an addi- 
tional factor iSi 1 '. The contribution to the bilinear action 
(~ V 2 ) is given by ^(Sy (<S^) 2 ) , to the quartic terms 



m 



TO + 1 



(m + l)(m + 2) 7A 3 



E 



to(to — 1) 
c k m{m + 1) 



4-fe k 

lk=o e - € + 



to(to — 1) (m + 1)(to + 2) 



+A 



k=l 
2 



e 4 ~ k e k . 



E 



dkm(m + 1) 
e 5 " fe d 



_fe=0 
(m + 1)(to + 2) 




to(to — 1) 

c 5-kh 



H 



(0) 



c 5-fc k 



9/x 



t 



2 <V 



(22) 



where cq — 3m/(m + 1), c\ = — C3 = 2, c 2 = 0, C4 = 
— 3(to + 1)/to, rf = 2c , d\ = d 4 = 3, d 2 = d 3 = —1, d 5 = 
— 2C4, e\ = 3, e 2 = e 3 = 4, and e_ 2 = (3 — 2m)U + 2/1, 
e + 2 = (1+2to)C/— 2JL. Finally, the replica mixing disorder 
vertex reduces to 
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In the absence of the speckle, G = 0, and the effective 
action (|10[) simply corresponds to n identical copies of 
the action for the clean system as given in section III Al 
and as derived previously.- Switching off the coupling of 
the disorder to the hopping matrix elements (7 = 0), we 
recover the case of pure uncorrelated potential disorder. 
In fact, for 7 = 0, the coefficients derived in this section 
are found to be identical with the ones derived in sec- 
tion [TIB] demonstrating that our results are independent 
of the order in which the Hubbard-Stratonovich transfor- 
mation and the replica disorder average are implemented. 
Whereas for pure on-site disorder, it is more economi- 
cal to start with the Hubbard-Stratonovich transforma- 
tion, the presence of simultaneous off-diagonal disorder 
enforces us to first average over the disorder and to per- 
form the more tedious strong-coupling expansion after- 
wards. 

Interestingly, even in the presence of simultaneous cor- 
related hopping disorder (7 > 0), the disorder vertex G 

vanishes as K± = determining the positions of the 
Mott-lobe tips in the clean system. At least to order A 3 , 
at such commensurate fillings, the only effect of disorder 
is to renormalize the coefficients of the clean system sug- 
gesting that a direct MI to SF transition obtains in the 
presence of weak disorder. 



III. RENORMALIZATION GROUP 

To investigate the effects of weak disorder on the phase 
diagram of the BH model and in particular to study the 
instability of the MI state towards the formation of a 
BG, we analyze the effective field theory ([14]) by means of 
the renormalization group (RG) method. In this section, 
the one-loop RG equations are derived. Further, it is 
explained how the different phases can be identified from 
the RG flow. 



A. Derivation of RG equations 

For convenience, we rescale time and length to dimen- 
sionless units as Xq — Ut and Xi — Ar^ (i — 1, . . . ,D), 
respectively, where A is the initial momentum cut-off of 
the theory. Note, that A is not sharply defined since we 
have performed a naive coarse graining in order to take 
the continuum limit and to invert the hopping matrix in 
the long- wavelength limit. In any case, one should ex- 
pect A = 2ir/a with a of order unity or a few lattice 
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FIG. 2: a) Diagrammatic representation of the interaction 
vertex h which is diagonal in the replica index and local 
in imaginary time and the replica mixing disorder vertex g, 
which depends on two different imaginary times, b) Diagrams 
contributing to the renormalization of the mass coefficient r. 
Note that the contraction of the disorder vertex also gives rise 
to a renormalization of the temporal gradient terms 71 and 
72 as described in the text. In c) and d) the diagrams leading 
to a renormalization of h and g, respectively, are shown. We 
restrict ourselves to one- loop order. Diagrams vanishing in 
the replica limit n — > are not shown. 



constants. In the dimensionless units, the momentum 
cut-off is A = 1. In addition, we rescale the fields as 
ip a = ilJ a \/K 3 A 2 - D /U, to obtain 



S c ff 

So 

Sa 



So+S a +S h , (24) 
/ (k 2 - ijioj + + r) \va{k,uj)\ 2 , 

/ dx o / d £ 'a;|^ a (a;,a;o)| 4 , 

J d^o J dx' J d c ' x\ip a (x, x a )\ 2 \ip b (x, x' )\ 2 , 



with dimensionless coupling constants 71 = K\U /(A 2 K 3 ) 
and 72 = K2U 2 / (h 2 K 3 ) for the temporal derivative 
terms, 73 = 1 for the spatial gradient term, and 
r = R/(A 2 K 3 ) for the mass coefficient. The rescaled 
interaction and disorder strengths are given by h = 
HU/(A 4 ~ D K 2 ) and g = G/(A 4 - D Kl). We wrote 
the bilinear actions Sq in frequency and momentum 
space using the Fourier transformation <p a (x,xo) = 
f, <p a (k, u})e~~ l ( kx+UT \ For brevity, we have defined 

k = i^r D /|*|<i d D fc and I = (27T)- 1 Ho d«. Note 
that since we focus on the zero-temperature limit, we can 
work with continuous frequency integrations rather than 
with discrete Matsubara summations. 

In order to derive the RG equations, we successively 
eliminate modes of highest energy with momenta on the 
infinitesimal shell e~ d/ < \k\ < 1 and consecutively 
rescale momenta, frequencies, and fields. To do so, we 
decompose the fields, 
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ip a (x 7 x ) = ip<(x,x ) + ip>{x,x ) (25) 
= f [ $ a {k,u)e~^ kx+ ^ 

J\k\<e~ dl Jui 

Je- dl <\k\<l Ju 

into 'fast' fields ip^ containing modes from the infinites- 
imal momentum shell and 'slow' fields tp£ depending on 
the remaining momenta. Decomposing the action ac- 
cordingly and taking the trace over the 'fast' fields (f> , 
we find that the corrections of order d^ to the action 
^If = S< + S< + S< are given by dS< = (S<> + 

5<>) ,>-i((5<>+5<>) 2 )^ where S<> andS<> con- 
tain all interaction and disorder vertices with two slow 
and two fast fields. The superscript on the average over 
the quadratic terms indicates that only connected dia- 
grams contribute. The diagrams giving rise to a renor- 
malization of order dl of the coefficients of the bilinear 
action, as well as the interaction and disorder vertices are 
shown schematically in Fig. [5] Using the correlator 



fields (ifia — ► ( Pae~ Xdl ), we obtain the renormalization- 
group equations for the coupling constants of the bilinear 
action, 

dv — 

— = 2r + 2I 1 h + I g, (29a) 
dl 

^ = -(4-2)7i-/ 2 7i5, (29b) 

^ = ~2(d z - 1) 72 -I 2 (/o7i 2 + 72)5, (29c) 

where we have set the scaling dimension of the fields to 
A = (D + d z + 2)/2 to ensure that the coefficient 73 = 1 
of the spatial gradient term is not renormalizcd. Further, 
we have defined h = 2S D /(2n) D h and g = 2S D /(2n) D g 
to absorb pre-factors. 

The quadratic corrections —i((iS^ > + iS< > ) 2 corre- 
spond to contractions of two vertices and lead to a renor- 
malization of h and g. The corresponding diagrams are 
shown in Fig. [2b and Fig. [2ji, respectively, and are easily 
calculated by using Eq. 1261 Including the contributions 
from rescaling the RG equations for the interaction h and 
the disorder strength g are given by 



( V l(k,u)<p b (k',u/)) = (27r) D+1 S ab S(k-k') 

x6(Lu-Lu')C (k,r), 
C {k,Lu) = (fc 2 -i 7l w + 72 w 2 + r)- 1 ,(26) 

we obtain immediately 

(5i<>)o ' >=4 (^ /lMl ^lf l ^ (fc ' w)|2 ' (2?) 

with Sd the surface of the D-dimcnsional unit sphere and 
the inner frequency integral I\ = ^Co(k — 1,oj) = 

I-co ^ C *o( k = = [4(1 + r) l2 + °7f]- 1 / a for 72 > 
and Ix = for 72 < 0. Comparing Eq. [57] with the bilin- 
ear part So of the effective action (|25|) . one immediately 
recognizes that the contraction of the interaction vertex 
leads to a renormalization of the mass coefficient r. Due 
to the presence of two independent times in the disor- 
der vertex, the inner frequency of the 'fast' fields in the 
diagram (S< > )o,> is not free and we obtain 

(<5 S <> K> = 2 7^I^E / rCo(l,^ a (fc,c)| 2 . 

\ Zn ) a Jul Jk 

Therefore, by expanding for small frequencies, Cq(1,lu) ~ 
h ~ Ifai(-iuj) - (7 2 72 + I 3 7i> 2 with I Q = 1/(1 + r), we 
find that the disorder vertex not only renormalizes the 
mass but also the coefficients 71 and 72 of the temporal 
gradient terms. 

From the first-order vertex corrections (5^ > +5^ ;> )o,> 
shown schematically in Fig. [3b>, and from the rescaling 
of frequency (ui — > u>e d - dl ), momenta (k — ► fce d/ ), and 



— = -{D + d z - A)h - 5I 2 h 2 - 6I$gh, (29d) 

^ = -{D-A)g-2{ll + h)g 2 ~4I 2 gh, (29e) 

where h := JZo ^ (1, w)C *(l, w) = I h/2 for 72 > 
and h = Zo/(2|-yi[) for 72 < 0. 

The above frequency integrals I\ and J 2 are defined for 
r+1 > which holds in the MI phases as well as in the SF 
sufficiently close to the insulating states. Further, the RG 
equations are valid only for sufficiently weak interaction 
h and disorder ~g. Note that since the effective field theory 
is dual to the initial Bose-Hubbard model, the smallness 
of the interaction vertex is guaranteed by t/U -C 1. 

B. Numerical Integration 

In the following, we integrate the RG equations [29k - 
d numerically whereas the initial values of the coupling 
constants 71 (0), 72(0), r(0), h(0), and g(0) are deter- 
mined by the microscopic parameters of the disordered 
BH model as explicitly derived in section [TTJ In the re- 
mainder of this section, we describe how the different 
phases of the disorder BH model are identified from the 
scale dependence of the coupling constants. The result- 
ing phase diagrams are presented in section HVl 

Since 71 (0) 7^ except for special values oifi/U mark- 
ing the tips of the Mott lobes, we use a dynamical expo- 
nent d z = 2. Although away from the tips, 72 eventually 
decreases exponentially under the RG, it is important to 
include 72 since the initial ratio 72 (0) /71 (0) becomes ar- 
bitrarily large in the vicinity of the tips, thereby strongly 
modifying the RG flow on small scales. 
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Whereas on the mean-field level, the phase boundary is 
determined by the sign change of the bare mass coefficient 
r(0) this parameter is renormalized by the interaction h 
and the effective disorder g leading to a shift of the phase 
boundary between the insulating and superfluid phases. 
In the insulating phases, r(l) — > oo, corresponding to a 
freezing of the system on large scales, whereas in the SF 
phase r(l) starts to diverge to negative infinite values. 
In the latter case, we stop the integration at a scale l s 
where r(l s ) + 1 = 0+ and the frequency integrals Ii be- 
come singular. Note that the singularity indicates that 
our approach is valid only in the strong coupling regime. 
As soon as the system flows to weak coupling, the RG be- 
comes unstable. However, sufficiently close to the local- 
ized states it is possible to estimate the superfluid density 
as pi ~ -r(l s )/h(l s ) exp(-Dl s ). 

The instability of the MI towards the formation of a 
BG is indicated by the divergence of the effective disor- 
der strength g. Since the RG equations are valid only in 
the weak disorder regime, we stop integration at a scale 
I* where — 1 yielding an estimate £ ~ aexp(7*) for 

the correlation length corresponding to the typical linear 
dimension of Mott insulating islands. Note that the scal- 
ing dimension of ~g indicates the relevance of the effective 
disorder in D < 4 and hence in any real physical system. 
Therefore, even for an infinitesimally small initial value 
g~(P), the system will eventually become unstable towards 
the formation of a BG although this might happen on as- 
tronomically large scales. In a finite system of linear di- 
mension L — <zexp(/ max ), the system looks ordered and is 
indistinguishable from an incompressible MI if / max < I*. 
In the following, we will investigate systems of different 
sizes, comparable to typical dimensions of optical lattices 
but also on the order of much larger condensed matter 
systems, and determine the crossover between the MI and 
the BG using the condition Z max = /*. 



IV. RESULTS 

In the following, we determine the zero-temperature 
phase diagrams of the disordered BH model on a cubic 
lattice {D = 3). For any set of parameters of the lattice 
model, chemical potential /U, hopping t, on-site repul- 
sion U, and disorder parameters A and 7, we obtain the 
coupling constants of the effective field theory from the 
relations derived in Sec. |TTJ Note that the effective field 
theory provides a valid description of the long-wavelength 
physics only in the regime of strong-coupling (t/U <C 1) 
and weak disorder (A/J7 <C I). By analyzing the RG 
flow of the effective coupling constants, we are able to 
identify the different phases delineated in Sec. IIIII 

As discussed in Sec. Ill A[ in the absence of disorder and 
at incommensurate fillings, the phase transitions from 
the MI states to the SF are of mean-field type. However, 
although the interaction vertex h is irrelevant and de- 
creases exponentially under the RG, it renormalizes the 
mass coefficient r to larger values as can be seen from Eq. 
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FIG. 3: (Color online) Upper left: Phase diagram for the clean 
system showing the m — 1, m = 2 MI, and SF phases. The 
SF density p B is shown as a color gradient and, for compari- 
son, the mean-field phase boundaries are indicated by dashed 
white lines. Other panels show the increase of the BG regions 
with the system size L = aexp(Z max ) for pure on-site disorder 
with A/U = 0.15 . Along the solid white lines, the effective 
disorder vanishes, G = 0. Grey areas are not accessible within 
the present strong-coupling approach. 



I29al leading to a stabilization of the MI. In the upper left 
panel of Fig. [3j we show the increase of the Mott lobes 
due to the renormalization effect. 

We continue with the case of pure on-site disorder 
(7 = 0). Already on a mean-field level, the presence 
of disorder leads to a destabilization of the Mott lobes 
since the bare value of the mass coefficient is reduced 
by disorder as can be seen from Eq. HU and K { J? ] > 0. 
In Fig. [31 the phase diagrams obtained from the nu- 
merical integration of the RG equations for a disorder 
strength A/U — 0.15 are shown as a function of fx/U 
and zt/U . In addition to the decrease of the Mott lobes 
due to the mass reduction, disorder leads to an instability 
of the MI state towards the formation of a BG. In Fig. 
[3j the finite-size crossover between the two states for sys- 
tems with different linear dimensions L = aexp(7 max ) is 
shown. Whereas for small systems comparable to typical 
optical lattice dimensions (7 max = 4), the BG emerges 
only in a thin sliver close to the SF, the BG regions in- 
crease with the system size and almost completely take 
over the MI regions for l max = 12. This behavior is symp- 
tomatic of the fact that the BG phase is determined 
by exponentially rare compressible regions.— While the 
BG regions increase with the system size, the different 
Mott regions become thinner and thinner and are cen- 
tered around (]J/U) m = \J m(m + 1) — 1 marking the 
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FIG. 4: (Color online) Enhancement and shift of the m = 1 
MI phase due to correlated hopping disorder (right) compared 
to the pure on-site disordered case (left). The BG region is 
slightly enhanced due to the disorder in the hopping. Solid 
white lines indicate the vanishing of the effective disorder, 
dashed lines the corresponding mean-field (MF) phase bound- 
aries determined by Rf-°'(n,t,U) = 0. For the BG phase 
boundary, no MF prediction is available. 
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tip positions of the Mott lobes in the clean system deter- 
mined by i^ "* = 0. Along these lines, the disorder vertex 

G = — ^(k[°') 2 A 2 vanishes. Consequently, the tips are 
impervious to disorder (at least at one-loop order) and 
a direct transition from the MI to the SF obtains in the 
thermodynamic limit. 

The effect of simultaneous correlated hopping disorder 
Stij — 7(e, — ej) 2 is illustrated in Fig. [4] for a speckle 
intensity of A/U — 0.4. Compared with pure on-site dis- 
order (7 = 0), already a small coupling zjU = 0.1 leads 
to a significant stabilization of the MI state as shown 
for the m = 1 lobe in the close vicinity of the tip. In- 
terestingly, the Mott lobes increase in size and simulta- 
neously shift to larger values of the chemical potential 
~fi/U . Further increase of 7 leads to further increase of 
the lobes whereas the effect is even stronger for larger 
fillings to. This behavior can be understood from in- 
specting the dependence of the mass coefficient R on the 
disorder parameters (Eq. [21"j) . Whereas potential dis- 
order always leads to a decrease of the bare mass and 
therefore to a shrinkage of the Mott lobes, simultane- 
ous correlated hopping disorder of the form studied here 
changes the bare mass coefficient by AK[°\A 3 /t. Recall 
that K[ ^ = —dR^/d/j, changes sign at the tip positions. 

For fillings below the tip of the Mott lobe, K± < 0, and 
therefore the Mott state is further destabilized. In con- 
trast, above the tip, k[ ^ > leading to an increase of the 
bare mass and hence to a stabilization of the Mott state 
as compared to pure on-site disorder. Although these 
arguments neglect renormalization effects and are based 
on the mean-field picture, the comparison with the phase 
diagrams obtained from numerical integration of the RG 
equations (see Fig. [4j shows that renormalization effects 
lead to an almost constant shift of the mean-field phase 
boundaries to larger values of zt/U, leaving intact the 
intuitive mean-field picture. 



FIG. 5: (Color online) Phase diagrams as a function of speckle 
intensity A/U and relative hopping strength zt/U for val- 
ues Jl/U — 0.3 (left panel) and Jl/U = 0.45 (right panel) of 
the chemical potential, corresponding to fillings below and 
slightly above the commensurate value (p/U)i, respectively. 
In the upper row, the phase diagrams for pure on-site disor- 
der are shown and contrasted with the ones in the presence 
of simultaneous correlated hopping disorder with coupling 
Az'yU — 1 (lower row). The crossover between the MI and 
the BG is shown for a system size corresponding to Z max = 8. 
Note that a reliable determination of the BG/SF transition is 
not possible within our approach. Alternative phase bound- 
aries consistent with recent optical lattice experiments^ are 
indicated as dashed lines. 

Note that the phase diagrams in Fig. [4] are ob- 
tained by numerical integration up to a logarithmic scale 
^max = 5 corresponding to relatively small systems of 
roughly 150 3 lattice sites. Whereas for pure on-site dis- 
order the crossover from the MI to the BG sets in only 
very close to the superfluid phase, away from commensu- 
rate fillings (see also Fig. [3]), the presence of simultaneous 
on-site and hopping disorder leads to a sizeable enhance- 
ment of the BG region as shown in Fig. [4] for z lU = 0.1. 
As in the case 7 = 0, the BG regions increase with the 
system size but never touch the lines at (jl/U) m corre- 
sponding to the Mott tips of the clean system. As can 
be seen from Eq. [23l the hopping-disorder contribution 
to G is proportional to and therefore, regardless of 
the values of A and 7, the disorder vertex vanishes on 
approaching commensurate fillings. 

To further investigate the evolution of the different 
phases as a function of the disorder strength and to fur- 
ther illustrate the effects of simultaneous correlated hop- 
ping disorder, we map out the phase diagrams as a func- 
tion of the relative hopping strength zt/U and the speckle 
intensity A/U for different values of the chemical poten- 
tial ~fi/U . In Fig. O such phase diagrams are shown for 
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Jl/U = 0.3, smaller than the value (jl/U)i = y/2—1 at the 
tip of the first Mott lobe and Jl/U = 0.45, slightly above 
the tip. For both values, the phase diagrams for pure on- 
site disorder are compared with the ones calculated in the 
presence of correlated hopping disorder with izjU = 1, 
where the system size has been fixed as Z max = 8. 

Since the values of the chemical potential correspond 
to incommensurate fillings, the bare disorder vertex 7?(0) 
is different from zero and consequently, on a certain scale, 
g(l) will eventually become of order unity signaling the 
instability of the MI towards a BG. Therefore, a direct 
MI-SF transition is impossible in the thermodynamic 
limit. In a finite system, the disorder strength A/U has 
to exceed a certain threshold to induce a BG separating 
the MI from the SF (see Fig. This critical disorder 
strength decreases with the system size but also depends 
on the value of k[ ^ and hence, on the value of Jl/U. 

On approaching the commensurate value (Jl/U)i, k[ ^ 
decreases leading to an increase of the critical disorder 
strength necessary to induce a BG. This behavior can be 
clearly seen in Fig. both for pure on-site disorder and 
in the presence of simultaneous hopping disorder. For 
Jl/U = 0.3 the critical disorder strength is much smaller 
than for Jl/U = 0.45 which is only slightly above the 
commensurate value. Note that in the latter case, for 
7 = the MI-SF remains direct up to the largest value 
A/U = 0.4 considered here. 

The effect of simultaneous hopping disorder is in agree- 
ment with the trend suggested by Fig. |4j For Jl/U = 0.3 
below the commensurate value (Jl/U)i the phase tran- 
sition to the SF is shifted to smaller values of zi/U, 
whereas for Jl/U = 0.45 > (JZ/[/)i the localized phases 
are stabilized by a finite 7. Moreover, it can be clearly 
seen that in the presence of correlated hopping disorder 
the system is more susceptible towards Bose-glass forma- 
tion. For all values zt/U in the MI regime the crossover 
to the BG sets in at smaller values A/U as compared to 
pure on-site disorder. 

Interestingly, starting from the MI at t/U close to the 
SF phase, weak disorder can induce a transition to a SF 
rather than to a BG. Note that such a crossover is only 
possible in a finite system. As expected, sufficiently deep 
in the MI phase disorder induces the crossover to a BG. 
On further increase of the disorder strength we find a 
transition from a BG to a SF in agreement with phase 
diagrams obtained by stochastic mean-field theory^ and 
Quantum Monte Carlo simulations in D = 3^ but in 
disagreement with recent optical lattice experiments^ 
which seem to indicate that a disorder induced transi- 
tion from an insulator to a SF is impossible. In fact, 
the transition between the BG and the SF can not be 
predicted within our approach. Since both the effective 
field theory and the RG analysis are valid only in regime 
of strong-coupling and weak disorder, we can only deter- 
mine the instabilities of the MI, either towards the SF 
controlled by a flow to weak coupling, or towards the 
BG controlled by a flow to strong disorder. Alternative 
BG/SF phase boundaries which would be consistent with 
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FIG. 6: (Color online) Evolution of the SF density as a func- 
tion of zt/U (left) and A/U (right) for pure on-site disorder 
(red) and with correlated hopping disorder (blue). 



the recent optical lattice experiments^ are indicated by 
dashed lines in Fig. O 

Finally, we analyze the dependence of the SF density 
p s on the effective hopping strength zt/U and the speckle 
intensity A/U for pure on-site disorder and in the pres- 
ence of simultaneous hopping disorder. For simplicity, 
we focus on the point (Jl/U)\ marking the tip of the first 
Mott lobe in the clean system. In the panel on the left in 
Fig. [6l we show the onset and increase of p s as a function 
of zt/U for different values of A/U. For larger values 
of A/U, the transition shifts to smaller values of zt/U 
signaling a shrinking of the insulating region. For all val- 
ues of A/U the correlation 7 leads to a reduction of p s 
being more pronounced for stronger speckle intensities. 
Whereas close to the transition, p s increases with A/U 
(Fig. O right-hand panel) we find p s to decrease suffi- 
ciently deep in the SF. Simultaneous, correlated hopping 
disorder, inherent to an optical speckle, enhances this 
effect and moves the point where p s starts to decrease 
closer to the transition. 



V. SUMMARY AND DISCUSSION 

To summarize, we have studied the weakly disordered 
BH model on a cubic lattice in the strong-coupling regime 
and at zero temperature. By combining a disorder aver- 
age using the replica trick with a strong coupling ex- 
pansion around the localized limit based on a Hubbard- 
Stratonovich transformation, we have derived the corre- 
sponding long-wavelength field theory, providing explicit 
expressions relating the effective coupling constants to 
the microscopic parameters of the disordered BH model. 
From the RG flow of the coupling constants, we have an- 
alyzed the instabilities of the MI state, either towards 
the formation of a SF characterized by a flow to weak 
coupling, or towards a localized BG phase as indicated 
by a flow to strong disorder. 

The techniques presented here can be applied to vari- 
ous forms of simultaneous on-site and hopping disorder. 
Also the inclusion of finite temperature effects or an RG 
analysis beyond one-loop order are in principle possible. 
In fact, the latter is necessary to capture the BG-SF tran- 
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sition at commensurate fillings. 

Whereas the presence of simultaneous correlated hop- 
ping disorder as induced by optical speckle potentials 
does not change the overall topology of the phase dia- 
grams, we have found that in comparison to pure on-site 
disorder (7 = 0), the insulating phases are enhanced con- 
siderably, wheras the SF is destabilized. Moreover, the 
system becomes more susceptible towards the formation 
of a BG. 

Our phase diagrams are in qualitative agreement 
with stochastic mean-field theory^ and Quantum Monte 
Carlo simulations in D — 3 in the weak disorder regime, 9 
although the latter seem to disagree on the possibility of 
a finite-size crossover between the MI and the SF in the 
presence of weak disorder. Unfortunately, at present, no 
complete experimental phase diagrams are available. In 
particular, it has not been possible yet to distinguish be- 
tween the two insulating phases, the BG and the MI^ 
In contradiction to our and previou o 9 ' 19 theoretical work 
the experiments seem to suggest that a disorder driven 
transition from an insulator (BG or MI) to a SF is not 
possible. However, a controlled calculation of the phase 
boundary between the BG and the SF is not possible 
within our present approach. The experimental observa- 
tion that sufficiently deep in the SF phase the superfluid 
density decreases with the speckle intensity 2 ^ is in qual- 
itative agreement with our predictions. 

One should bare in mind that a direct comparison of 
our theoretical phase diagrams with the data obtained 
from optical lattice experiments suffers from various po- 
tential problems. First of all, the typical dimension of 
optical lattice systems is much smaller than of realistic 
condensed matter system and therefore one should expect 
much broader crossovers between different phases. A sec- 
ond complication is the presence of an optical trap po- 
tential leading to inhomogencous densities which might 
cause phase separation in the system or destroy commen- 
suration effects which are crucial for understanding the 
BH model. Third, the experiments suffer from heating by 
the lasers leading to temperatures which are not neces- 
sarily small compared to other experimental parameters 
as the effective Hubbard repulsion U . Also temperature 
changes during the experimental preparation and mea- 
suring processes are likely to obscure the data. 

However, despite all these challanges, the optical lat- 
tice experiments have been billed as the ultimate quan- 
tum simulators of paradigmatic lattice Hamiltonians such 
as the one defining the disorder BH model. In this work, 
we have presented thermodynamic quantum phase dia- 
grams for this model, where we have used disorder distri- 
butions and correlations as induced by an optical speckle 
lens to arrive at a comprehensive picture. 

Our major finding that the disorder vertex vanishes to 
leading order at commensurate boson fillings provides a 
plausible explanation for the long-standing controversy 
of wether a direct MI-to-SF is possible in the presence 
of weak disorder. On approaching the tips of the Mott 
lobes, the bare value of the disorder vertex is strongly 



suppressed and consequently the scale on which the BG 
is observed becomes enormously large. Not surprisingly, 
in earlier Quantum Monte Carlo (QMC)i&±ii±ii±&±£ and 
density matrix renormalization group studies^ of rel- 
atively small systems, no indications of a BG have 
been found at or close to commensurate fillings whereas 
from more recent QMC data on larger two£ and three 
dimensional 9 systems it has been concluded that even at 
the tips, a small amount of disorder is sufficient to induce 
a BG disrupting the MI-SF transition. However, based 
on numerical simulations it is impossible to answer the 
question if in the thermodynamic limit a direct MI-SF 
transition obtains in the presence of infinitesimally weak 
disorder. 

Various analytical methods^^^^ ' 12 ' 15 i 19 have been 
employed to find an answer to this question, arriving at 
different conclusions. While there is consensus on the 
absence of a direct MI-SF transition at incommensurate 
fillings, with the exception of the stochastic mean-field 
theory 19 which presumably gives wrong results below the 
upper critical dimension D = 4, the case of commensu- 
rate boson fillings remains highly controversial. In their 
seminal paper, Fisher, et ali, who raised the question 
in the first place, argued that a direct MI-SF transition 
at commensurate fillings was unlikely, though not funda- 
mentally impossible. A work based on strong-coupling 
series expansions^ arrived at the same conclusion based 
on the observation of Lifshitz rare regions. In the one- 
dimensional case, Svistunov^ demonstrated by an RG ar- 
gument that commensuration is always asymptotically 
destroyed by disorder and hence that the transitions 
at commensurate and incommensurate fillings are con- 
trolled by the same BG fix-point. Based on a double 
epsilon expansion, Weichman and Mukhopadhyay 2 ^ sug- 
gested the existence of a universal disordered BG fixed- 
point in higher dimensions, irrespective of the filling, at 
which particle-hole symmetry is statistically restored. Fi- 
nally, in the limit of large fillings, Herbut^^ demonstrated 
the absence of a MI-SF transition at commensurate fill- 
ings by using duality mappings to the sine-gordon model 
in D = 1 and the three dimensional Higgs electrody- 
namics in D = 2. Whereas all these works support the 
destruction of the direct MI-SF by infinitesimally weak 
disorder, using real-space RG technique o 12 ' 15 a direct MI- 
SF transition at commensurate fillings up to a critical 
disorder strength is found. 

In our present work, we find that at incommensurate 
fillings already an infinitesimal amount of disorder is suf- 
ficient to induce a BG separating the MI from the SF, 
in agreement with previous results^ 2 -'^£&£ ' 12 i 15 . The 
vanishing of the disorder vertex g at the tips of the Mott 
lobes to leading order (A/U) 2 clearly indicates that com- 
mensuration effects are crucial for understanding how the 
system couples to an infinitesimal amount of disorder. 

Whereas on a one-loop level, consistently capturing 
disorder contributions of relative order (A/U) 2 , a di- 
rect MI-SF transition obtains at commensurate fillings, 
this result is unlikely to hold at the level of a two-loops. 
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First, the next order in the disorder cumulant expan- 
sion yields a contribution of order (A/U) 4 to the disorder 
vertex g which does not vanish at the tips. A contribu- 
tion of the same order comes from the renormalization 
by the term ~ A 2 [(p* a (T)d T (p a {T)][(pl(T')d T (p b {T')] at one 
loop order ^ which is less relevant than the generic dis- 
order vertex but has to be included at the tips as the 
leading order quartic disorder vertex. Such a term arises 
in the replica theory from the coupling of potential dis- 
order to the linear time derivative term. Therefore, g 
does not remain zero at the tips and the MI-SF is not 
protected by symmetry. However, including this term in 
a one-loop calculation is no longer consistent as a term 
~ A 2 |(^ a (T)| 4 |<^6(r')| 4 , generated by the coupling of dis- 



order to the interaction vertex, gives rise to a renormal- 
ization of g of order (A/U) 2 at two loop order. Con- 
sequently, a two-loop calculation seems unavoidable to 
settle this issue. We believe that an extension of the cal- 
culations presented here in combination with the double 
cpsilon expansion by Wcichman and Mukhopadhyay 28 
allowing to identify potential disorder fixed points will 
provide conclusive answers, not only to the existence of 
the BG but also on the nature of the transitions. 
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